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SUMMARY 


Computer-generated graphics in real-time helicopter simulation produces objectionable scene- 
presentation time delays. In the flight simulation laboratory at Ames Research Center , it has been 
determined that these delays have an adverse influence on pilot performance during aggressive tasks 
such as nap-of-the-Earth (NOE) maneuvers. 

Using contemporary equipment , computer generated image (CGI) time delays are an unavoidable 
consequence of the operations required for scene generation. However , providing that magnitude 
distortions at higher frequencies are tolerable , delay compensation is possible over a restricted fre- 
quency range. This range , assumed to have an upper limit of perhaps 10 or 15 rad/ see, conforms 
approximately to the bandwidth associated with helicopter handling qualities research. 

A compensation algorithm is introduced here and evaluated in terms of tradeoffs in frequency 
responses. The algorithm has a discrete basis and accommodates both a large , constant transport 
delay interval and a periodic delay interval , as associated with asynchronous operations. 


INTRODUCTION 


Flight simulation research at the Ames facility includes 
various computer, motion base, and subsystem configuration 
options. Either an electro-mechanical camera/model board 
device, or a Singer-Link CGI system are available for scene 
presentation. Despite an adequate bandwidth capacity for 
helicopter research, the model-board's display leaves some- 
thing to be desired because of its restricted field of view. In 
contrast, the multi-window CGI, while manifesting minor 
deficiencies in scene content and texture, is valuable for pilot 
cues. However, owing to delays in the transmission of posi- 
tion and orientation information, CGI frequency responses 
are degraded within the operational bandwidth of the pilot. 

The Singer-Link CGI system operates asynchronously 
with a host computer, although the “time of data transfer” 
is tagged by an accompanying interrupt signal. Data flows 
sequentially through three CGI processors: (1) the Perkin- 
Elmer Computer, (2) the Frame Calculator, and (3) the Scan 
Line Processor. This flow is outlined in figure 1. The pro- 
cessors have the following functions: 

Perkin-Elmer Computer 

The P-E computer is also called the scene management 
computer. One of its functions includes the collecting of 
mathematical descriptions of objects to be used in the scene 
which requires extensive data retrieval from disc storage. 
The arrival of new position and orientation data in the P-E 


computer is accompanied by an interrupt signal; however, 
scene management proceeds with previous data at a 30 Hz 
rate. 


Frame Calculator 

The Frame Calculator is also called the geometric proces- 
sing computer, or digital image generator (DIG). Geometric 
processing is the conversion of mathematical descriptions of 
three-dimensional data base objects into two dimensions, as 
associated with the coordinates at the display. 


Scan Line Processor 

The Scan Line Processor is also called the video processor. 
It performs the steps necessary to define the image at each 
picture element (including considerations of occlusion), 
and generates the video display signals. 

When these processors perform their required functions, a 
time delay occurs between the commanded presentation of 
the host computer and the actual display. This delay consists 
of “a pure transport delay” component plus a “periodic” 
component, which are quantified herein. The discrete nature 
of these delay sources suggests that a compensation scheme 
also have a discrete (rather than continuous) basis. This 
is the approach used here. Whenever possible, the results are 
displayed using conventional parametric representations, 
such as magnitude and phase diagrams. 



THROUGHPUT 

The major component of time delay is the workload 
required to perform the above processes. This workload is 
equivalent to the CGI system throughput: 

. the system throughput (defined as the 
time it takes to complete a scene after the 
parameters of the viewpoint become available) 
must be so fast that the delay does not dis- 
tract the observer. The commonly used 
update requirement is 100 msec” (ref. 1). 

“Throughput” is pure transport delay. However, other 
components of time delay also exist in simulation; the sum 
of these components degrades simulation fidelity. 

Delay Sources 

Five sources of time delay in discrete simulation are 
identified in this section. Two of these, pure transport delay 
and periodic delay, are the major causes of CGI response 
degradation. 

Time delay due to discrete realization- When the discrete 
sample rate is high compared to the frequency content 
delivered by the model, phase differences between con- 
tinuum and discrete representations are generally small. 
Equivalent time delays are functions of frequency, and over 
some bandwidth may actually produce lead rather than lag, 
depending upon the discrete parameters. This influence is 
relatively negligible for a “sufficiently band-limited” system 
(ref. 2), such as that assumed here. 

Loiter delay- In the context of this paper, loiter delay 
does not exist. Quoted values for loiter delay, on the order 
of 20 msec, indicate that it is usually identified with the 
mean value of the P-E computer’s asynchronous operational 
delay, more properly given as 16.67 msec. However, this 
quantity is not significant here; the host computer’s I/O 
cycle time T actually determines the frequency distribution 
function of the delay, and it is periodic rather than random. 
Periodic delay is handled below, and includes this “loiter” 
contribution. 

There is some evidence that a small transmission delay 
(which might be called loiter delay) of approximately 6 msec 
exists in our facility configuration (ref. 3), but this value is 
so small that it is ignored here. 

Pure transport delay - The P-E Processor, the Frame 
Calculator, and the Scan Line Processor operate in parallel at 
30 Hz. However, information flows serially through each of 
these components because of their pipeline architecture. 
Scene presentation thus occurs 100 msec after information 
utilization (not acceptance) by the P-E computer. This delay 


value is constant, and is given here by T c \ it is the major 
source of CGI response degradation. 

Periodic delay- Asynchronous operation of a CGI proces- 
sor with respect to a host computer produces a time delay 
that appears to be random, but is actually periodic. Harmon- 
ics of this delay may also be present, depending upon the 
host computer cycle time T. Periodic delay, given here by 
7^(m), varies at a 30-Hz rate. (This will be discussed under 
Maximum Time Delay.) 

ZerO’Order hold delay- Because of the 30-Hz data refresh 
rate of CGI monitors, scenes are presented in accordance with 
the zero-order hold phenomenon, with an average equivalent 
time delay of 16.67 msec. When this phenomenon is con- 
sidered along with the electronics of continuous scan line 
generation by the video processor, which also has an average 
time delay of 16.67 msec, the net result is a pure transport 
delay of 33.33 msec. The zero-order hold delay may thus be 
ignored when the total value of 33.33 msec is attributed to 
the Scan Line Processor. 

SCENE PRESENTATION DELAY 

Figure 1 outlines the relevant processes in the CGI loop 
for the purpose of examining scene presentation delays 
which are caused by pure transport delay and periodic 
delay. The host computer produces data each T, which is 
assumed to be correctly synchronized with real time, and 
immediately transferred (neglecting any residual loiter delay) 
to the P-E buffer (X). As shown in figure 1, this transfer is 
accompanied by an interrupt signal. This Data-Ready Inter- 
rupt (DRI) is transmitted by the host computer each I/O 
cycle time T, and accompanied by the data buffer X. The 
interrupt is serviced by the P-E computer without interfering 
with current operations on the extant buffer Y. The data 
buffer X replaces the data buffer Y every 33.33 msec within 
the P-E computer. 

Sample Rate Conversion 

Determining the maximum value of periodic time delay 
for an arbitrary host-computer cycle time T requires a 
numerical analysis based upon least common multiples. For 
cycle times restricted to integer multiples of milliseconds, the 
analysis proceeds as follows: 

The host computer produces a signal x(n) sampled at the 
rate F= 1/T , and the P-E computer observes the signal 
y(m) sampled at the rate F = 1/T\ Where N and M are 
integers, the ratio of sampling periods may be expressed as a 
rational fraction: 

T/T = F/F - M/N (1) 
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The ratios of integers (M/N) will be seen to be important in 
the following material. Using the P-E’s sample rate of 30 Hz, 
equation (1) becomes 

(33.33. ..)/!T = 100/(3 T) = M/N (2) 

where T is an integer number of milliseconds. Given that N 
is the least common divisor in equation (1), it may be shown 
(ref. 4) that the system response is periodic in m with period 
N. The period of time delay is thus given by 

P ~ NT' = 7V/30 (sec) (3) 

when the fraction given by equation (2) has been reduced to 
its lowest terms, i.e. M and N are relatively prime (to each 
other). Values for M and N over a wide range of host- 
computer cycle times T are shown in tables 1 and 2. 
For example, when T = 60 msec equation (2) reduces to 
M/N = 5/9. In the range of host computer cycle times from 
one to 70 msec the period P may be as small as 0.1 sec (e.g., 
when T = 50 msec), or as large as 100 times the host com- 
puter’s cycle time (e.g., 6.9 sec for T- 69 msec). The period 
of time delay is given in figure 2(a), where the cycle time T 
of the host computer is an integer number of milliseconds. 
Harmonics of this period also occur for various cycle times. 


Maximum Time Delay 

Both the maximum and average scene presentation 
time delays are shown in figure 2(b), as a function of cycle 
time. This figure illustrates the discontinuous behavior 
of the asynchronous system as a function of its primary 
parameter, T. 

The data in figure 2(b) agree with measured data acquired 
in a previous study (ref. 5), but this analysis draws a slightly 
modified conclusion. Paulk reports a linear relationship 
between T and lag, but it is shown here that only to a first 
approximation does a linear relationship exist between the 
visual time delay and the cycle time. Also, the additional 
computer cycle of the reference (a technique of measure- 
ment) is not pertinent. Without compensation of any sort, 
the model/CGI average time delay may be approximated by 
100 msec plus one -half of the cycle time of the host com- 
puter; but for 50 msec, for instance, this is a very pessimistic 
approximation. It is assumed that the host computer will 
transmit the correct command as synchronized with real 
time; otherwise the results are distorted by the discrete inte- 
gration algorithms. This assumption is not invalidated if the 
drive signal is a control input of the pilot. Inputs, which are 
generally force or moment proportional, may be assumed to 
be correctly integrated within the host computer. 

The periodic component of time delay may be written 
analytically. At the 30-Hz rate, where “time” may be given 
as t ~ m/ 30 sec (within the P-E computer), the periodic delay 


component is given by 

T A (m) = T A (m + U)= T{m/(30T) - lm/(30r)|f (4) 

In this expression the half brackets Q J) denote the least 
integer operation. The interval of time T A (m) constitutes 
the “projection interval” over which the P-E computer 
calculates extrapolated positions by use of the products of 
velocity times T A (m). This projection operation will be 
discussed further once the subject of compensation is 
introduced. (See “Projection.”) 

The total time delay in the asynchronous system of 
figure 1 has the periodic component given in equation (4). 
The minimum time delay is the pure transport delay 
T c = 100 msec, while the maximum time delay is given by 
T c + T^{m), where m increments at 30-Hz. Time delays as 
functions of cycle time are illustrated in figures 3(a) through 
3(1) for cycle times in the range of 49 to 60 msec (each 
msec). The fundamental period of time delay is also shown in 
each of these figures, and agrees with values previously given 
in figure 2(a) for a large range of cycle times. Harmonics are 
also observed for certain cycle times. The maximum value for 
T^(m) is less than the cycle time T, and the average value 
is one-half the maximum value. These facts are demonstrated 
in the Appendix. Tables 1 and 2 are provided for clarifica- 
tion. In these tables the value m(max) is the value of m 
when T A (m) achieves its maximum value over the period P. 

From the above it is concluded that total time delay is 
always less than the sum of cycle time plus pure transport 
delay. 

MODEL TRANSFER FUNCTION 

In this section a second-order transfer function is devel- 
oped that approximates specific aircraft behavior; in simu- 
lation this behavior was recently used to identify delay 
problems in CGI responses. This transfer function is selected 
as a basis for the study of the asynchronous logic of figure 1 . 

Scene presentation delays have been identified as causing 
degradation related to erroneous pilot cues in the XV-15 
(Tilt Rotor, UH-60A) simulation model, (ref. 6) especially 
in roll response (ref. 7). For the purpose of simplification, 
an approximate second-order system is used in this analysis. 
The system closely matches the observed hover response 
of the XV-15 in the roll axis. The second -order system 
has the natural frequency co = 7.259 rad/sec and damping 
£ = 0.4547. Using these parameters, the continuous response 
of the model is shown as a parametric curve in both 
figures 4(a) and 4(b). The discrete equivalent of this system 
is also given at the specific cycle time of 60 msec. 

In order to convert the system into an equivalent discrete 
model, the specific integration algorithms (ref. 8) used at 
Ames Research Center have been used, and are given here in 
z-transform notation. These algorithms consist of the Adams’ 
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second-order integration formula for the acceleration-to- 
velocity integration (a predictor), 


A CO = 


(T/2) (3z - 1) 
z(z - 1) 


( 5 ) 


and the trapezoidal formula for the velocity-to-position 
integration (a corrector). 


h CO = 


(T/2)(z + 1) 

Z - 1 


( 6 ) 


By using these algorithms and defining the parameters, 

A= coT 
a 0 = (4 S-A)A 
a, = 2A 1 - 1 6£A + 4 

a 7 = 3A 2 + 12M - 8 (7) 

the discrete system response function (z transform) is 
produced: 


F(z) = 


i4 2 (3z 2 + 2z - 1) 

4 z 3 + a 2 z 2 + 0 ^ +a 0 


( 8 ) 


Magnitude and phase responses of the discrete system given 
by equation 8 are a function of the cycle time of the host 
computer 7. The resultant curves are tagged “discrete simu- 
lation” in figures 4(a) and 4(b), where the cycle time is 
assumed to be 60 msec; they show that this relatively large 
cycle time is sufficient (small enough) for the approximation 
of continuum results for this transfer function out to about 
20 rad/ sec. 

To approximate the influence of CGI time delay on the 
system F(z), the total phase angles owing to the minimum, 
average and maximum delays are shown in figure 4(b) as 
computed by simple algebraic operations. (Time delays do 
not influence the magnitude responses of figure 4(a), but do 
produce large phase errors, as shown in figure 4(b)). 

Figures 4(a) and (b) have been obtained by using z trans- 
forms. They constitute a basis for comparison with other 
data in this document which has been obtained by use 
of time-series data. 


COMPENSATION 


In this section, modified position and velocity terms are 
developed that project positions (and orientations) much 
further than previous formulations. The assumptions of this 
projection are such that the entire delay interval is 


accommodated, with good frequency response characteristics 
over the entire operational bandwidth. Acceleration terms 
have been avoided in the development of this algorithm 
because: 1) they are not concurrent with velocity and 
position terms in real-time simulation, and 2) their high- 
frequency components, prior to “airframe attenuation,” are 
not generally consistent with the assumed band-limited 
behavior. 

In the host computer, operations by the integration 
algorithms (5) and (6) produce the vehicle velocity v(n) 
and position u(n ), which are transmitted to the X buffer of 
the P-E computer with an accompanying data-ready interrupt 
signal, as indicated in figure 1. Heretofore, only compensa- 
tion for the periodic interval T^im) has been used because 
the projection of an interval greater than this value cannot be 
supported under a purely linear hypothesis. It is known, 
however, that the additional pure transport delay component 
T c of 100 msec also exists prior to scene presentation. 
Indeed, for less exacting tasks than helicopter research, 
this delay seems to be generally accepted (ref. 1). Failure 
to account for this larger delay interval is due to the fact 
that general delay compensation does not exist. The com- 
pensation problem is put in its proper context in the “linear 
approach” section. 

Linear Approach 

“Projection” has been defined as the linear extrapolation 
over T^{m) that is performed within the P-E computer. 
The question arises as to the influence of also extrapolating 
positions over the pure transport delay interval within the 
host computer. The resultant pure linear compensation for 
the entire delay interval is here examined in terms of its 
component parts. 

The technique of linearly extrapolating position over the 
additional pure transport delay interval of 100 msec success- 
fully eliminates the phase angle error given in figure 4(b), 
but deteriorates the magnitude response of figure 4(a). 
These phenomena are illustrated in figures 5(a) and 5(b). 
where the discrete curves were obtained from a statistical 
analysis of time series data with Gaussian inputs. The time 
series data was created from a model of the asynchronous 
logic as shown in figure 1 . The use of the unusual cycle time 
7= 66.66 msec is related to the statistical analysis of (what 
appears to be) a nonstationary system, under the conditions 
of minimal data acquisition. 

Figures 5(a) and (b) each contain four parametric curves in 
order to illustrate the components of the extrapolation pro- 
cess, and the standard “projection” process. The continuous 
system curves are identical to those in figures 4(a) and 4(b). 
The curves labled “discrete, without linear compensation” 
represent the scene presentation without any attempt at either 
extrapolation or projection. Since this particular system does 
not use the standard, linear projection operation of the P-E 
computer, these curves approximate the average time delay. 
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As will be discussed under “Model Verification,” this average 
delay is 116.66 msec when T= 66.66 msec. The remaining 
two parametric curves in figures 5(a) and (b) show both the 
influence of linear extrapolation, or compensation, and the 
combined influence of linear compensation plus projection. 
The combined operations thus constitute the direct applica- 
tion of linear terms to the prediction of positional data. 

Figure 5(a) illustrates the basic problem encountered 
when linear compensation is used. As an extension of this 
approach, the lead-lag technique implies that sufficient lag 
must be superimposed to attenuate the magnitude error 
(linearly compensated, with projection) of figure 5(a), while 
retaining the phase characteristics of figure 5(b). Using 
linear analysis, recent work has been performed at Ames to 
“distribute the attendant system gain distortion (over fre- 
quency) so as to minimize gain -distortion effects on system 
responsiveness (co c ), pilot workload, and tracking accuracy 
(ref. 9).” This approach addresses the pilot -aircraft dynamics 
in the region of “crossover.” The nonlinear, double-tuned 
algorithm developed here is not an extension of this tech- 
nique, although the applicable bandwidth is assumed to 
contain the crossover frequency. 

Modified commands are developed below that account for 
the total delay interval, under assumptions that include 
tuning increments to both a specified frequency, and to the 
constant velocity condition. The resultant algorithm will be 
seen to retain its features in the interval between these 
discrete assumption points, and thus enhance the quality of 
position information (both translational and rotational) 
over the operational bandwidth. 


The Tuned Corrector 


by substituting u % = exp (jWkT c ). When equation (9) is 
tuned to a specific “cutoff frequency” W , the following two 
relationships are produced, 

(f 0 -f 2 ) tin U=(\ -cos U)IU (11) 

( f 0 + / 2 ) cos u + /i = sin u / u (12) 

where the cutoff bandwidth is given by 

U=WT C (13) 

The specific frequency W is defined as the upper limit of 
some “operational bandwidth.” The lower limit of this 
bandwidth is the zero frequency, or constant velocity condi- 
tion. By also tuning equation (9) to this condition, the 
additional relationship 


fo +/i +/ 2 = 1 (14) 

is produced. The simultaneous solution of equations (11), 
(12) and (14) then produces the intermediate coefficients 


fo 


U sin U - 2 cos U( 1 - cos U) 
~ 2 U sin U(l - cos U) 


(15) 


2 sin C/(sin U - U cos U) 
2 U sin U( 1 -- cos U) 


(16) 


In this section, a position increment is derived that 
accounts for the interval of pure transport delay. A sequence 
of velocity values is assumed available at 100-msec intervals. 
The compensation scheme is derived from this assumption, 
plus a selected parametric relationship. In the next section, 
a technique for approximating the required velocity sequence 
is given. 

Consider the three-parameter, velocity-to-position cor- 
rector relationship given by 


UsinU - 2(1 - cos U) 
2 U sin U( 1 - cos If) 


Using these coefficients, equation (9) represents the position 
transition from the point k to the point k + 1 . By forcing 
the point k to be coincident with the most recent data 
observation point of the host computer n , we may write 
equation (9) as, 


u k+ i =u k + (fo v k+i + fi v k + fi v k - 1 ) T c ( 9 ) 

In this expression time t is assumed to be incremented 
according to 

t = kT c (k = 0,1,2,...) (10) 

where T c is the large transport delay interval of 0.1 sec. In 
the following material, the “tuned” condition simply means 
that the stated condition satisfies the assumed relationship 
(eq. 9). For example, sinusoidal behavior may be described 


u k+i =u n + v ' T c < 18 ) 

where the applicable velocity v' for this increment has been 
given in terms of velocity values that are equally spaced^ by 
100-msec intervals, i.e., 

v= fo v k+i + fi v k + fi v k-i (19) 

A form of this corrector has also been developed using 
geometric relationships in an application to a system with a 
known dominant mode. It has been successfully applied to 
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the integration of aircraft states (ref. 10). In order to apply 
the corrector to an asynchronous system, however, velocity 
values symmetrically located about the point k (or n) with 
interval T must be developed. 


Symmetric Data 

The point k is defined as the most recent point of data 
delivery of the host computer. The point k + 1 thus defines 
a point 100 msec in the future, and the point k - 1 defines a 
point 100 msec in the past. Velocity values at these two sym- 
metric data points are estimated in this section using the 
actual data sequence within the host computer. 

In order to create symmetrical data points consistent with 
the above assumptions, sinusoidal prediction is used. Con- 
sider the velocity sequence {v w _ 2 , , v^j-of the host 

computer, where v n is the most recent value. In contrast to 
equation (10), time t is related to the host computer’s cycle 
time 7* by 

t = nT (n = 0,1,2,...) (20) 


Using these intermediate coefficients, the 
sequence for equation (19) is given by 

KTjfc) = v k + 1 =£o + ft sin U+g 2 cos U 

*0) = v k = v„ 

H-T/c) = vjc-i = So~gi sin U+g 2 cos U 
where U was defined in equation (13). 

The Algorithm 

In view of the above developments, equation (9) may be 
converted to a function of the velocity sequence of the host 
computer, 

u k+ i = u n **" v n-\ + 2 (30) 

where applicable coefficients for this difference equation are 
given by 


symmetric 

(27) 

(28) 
(29) 


A sinusoidal curve fit to the velocity sequence is given in 
the temporal space of the host computer by 


v(t - 2T) = go + g, sin [ W(t - 27)] + g 2 cos [ W(t -2 T)] 

( 21 ) 


where, for convenience 


v n -2 = K-27’) 


»«- i=v(-7) 


( 22 ) 


U + sin U( 1 - 2 cos V) (1 - cos U){ 1 + 2 cos V) 


2 W( 1 - cos V) 

+ 2 W sin V 

(31) 

(sin U - U)cos V 

(1 - cos U)( 1 + cos V) 


W{\ - cos V) 

W sin V 

(32) 

U - sin U 

1 - cos U 


| 

2 W( 1 - cos V) 

2 W sin V 

(33) 


v n = v(0) 

By defining the applicable bandwidth 

V= WT (23) 

the curve-fit coefficients may then be given by 


In addition to estimating position, the velocity at k + 1 
must be estimated for P-E projection purposes. This may be 
accomplished by using equations (24) through (27). The 
velocity that accounts for pure time delay (for later pro- 
jection purposes) may then also be written in terms of the 
host-computer velocity history, 




- 2v n _ x cos V + v 
2(1 - cos V) 


n -2 


(24) 


v k + 1 =doV n +d i v n . l +d 2 v n _ 2 
where the appropriate coefficients are given by 


(34) 


(1+2 cos V)v n - 2(1 + cos V)v n 1 + v n 2 

*' ISTf (25) 


ft = 


(1 -2 cos V)v n + 2 cos V-v n _ 2 
2(1 - cos V) 


(26) 


1 + (1 - 2 cos K)cos U i (1 + 2 cos F)sin U 

2(1 - cos V) + 2sin V (35) 

cos F(cos U - 1) (1 + cos F)sin U 

1 - cos V sin V (36) 
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1 - cos U sin U 

d 7 = + (37) 

2 2(1 - cos V) 2 sin V 

Since the b ^ and d ^ coefficients are constant for a given 
cycle time, the applicable compensation equations are the 
simple linear combinations given in equations (30) and (34). 
They constitute the real time workload of the host computer 
for each of the position and orientation commands. The 
workload of the P-E computer is not changed. It consists of 
extrapolation over the periodic interval, i.e., 

u ' = u k + 1 + T A( m >k + 1 ( 38 ) 

which has the same form as the original P-E algorithm. Note, 
however, that the symmetrical distribution of time delay 
about T c + 7^(avg) suggests that the compensation interval 
of the algorithm may be further extended to include the 
average value of periodic delay, with appropriate modifica- 
tion to the projection interval in the software of the P-E 
computer. Neither this extension nor the technique of 
accommodating the entire compensation algorithm within 
the P-E computer are discussed here. 

The algorithm is outlined in figures 6(a) and (b). The 
mainframe operations are shown in figure 6(a), and the 
P-E operations are shown in figure 6(b). The previously 
established second-order system is used in figure 6(a) 
for illustration, along with the standard integration 
algorithms (5) and (6). In figure 6(b), the “up arrows” refer 
to sample rate expanders (ref. 4), an operation also referred 
to as “interpolation.” If the host computer cycle time was 
less than 33.33 msec, these arrows should properly be 
pointing down in an operation called “decimation.” The 
“expanded” operations are asynchronous with those of 
figure 6(a), and produce periodic behavior in the quantity 
T a (z>z’). The z -transform operations referred to by z' have 
a 30-Hz update rate; the resultant time delay represented by 
( z ')” 3 is 100 msec. 

ASYNCHRONOUS MODEL 

For the purpose of computing frequency responses the 
pertinent asynchronous data flow of the CGI system was 
modeled. Interrupt logic was emulated by using time incre- 
ments with the granularity of 1/3 msec. Pertinent mainframe 
and CGI operations are as outlined in figures 6(a) and (b), 
where the z and z' transformations refer to different sample 
intervals, with periodic interaction. 

For statistical analyses, a worst-case cycle time of 
66.66 msec was selected for two reasons: First, this value 
approximates the value of 60 msec that was used in recent 
simulation studies, and is therefore currently relevant. 


Second, serious resolution bias errors (ref. 11) are avoided 
by using a sample rate that is a multiple of the P-E computer 
speed. Otherwise, the system appears to be nonstationary 
owing to its extremely long period, and conventional spectral 
analysis techniques yield dubious estimates of frequency 
response functions. 

Model Verification 

Figures 4(a) and (b) were obtained by using z-transforms 
as convoluted to frequency space, and thus represent a 
standard for verification of the asynchronous model, which is 
measured only by use of statistical analysis procedures. 

The estimated frequency response functions of figures 7(a) 
and (b) were obtained from the asynchronous model by 
using correlation and spectral analysis. Continuum results 
are also displayed in these figures, and parametric curves 
are presented from the analysis. In figures 7(a) and (b) the 
compensation algorithm was not used, so that the P-E 
computer model received the primary values given by u(z) 
and v(z) as shown in figure 6 (a). Two parametric curves are 
shown, indicating whether the clock in figure 6(b) was used 
(with projection) or not (without projection). 

The use of cycle time that is a multiple (two) of the 
cycle time of the P-E computer has minimized distortion 
in the spectral analysis so that the curves closely approximate 
the results previously given in figures 4(a) and (b). Using this 
particular cycle time, phase errors do not vary much from 
the continuum results owing to the fact that N = 2 in equa- 
tion (2) which results in a maximum T A (m) = 33.33 msec 
in equation (4). The average projection interval is only 
16.66 msec when T = 66.66 msec; hence, the phase angle 
history in figure 7(b) remains close to the “minimum bound- 
dary” given in figure 4(b). 

Projection 

The “projection” operation of the P-E computer does not 
appear to enhance performance in figures 7(a) and (b) to 
any great extent, but this is a misleading observation when 
viewed in frequency space. The small decrease in phase error 
is actually quite significant, and the resultant “smoothing” 
is quite important. Illustrations of this are made by con- 
verting to the time domain and using an accelerated sinu- 
soidal drive signal. This particular drive signal illustrates the 
bandwidth of interest because its frequency ramps from 0 to 
20 rad/sec in the displayed interval of 7 sec. 

To demonstrate the smoothing phenomenon in temporal 
space, the traces “without projection” and “with projection” 
are given in figures 8(a) and (b), respectively, using the 
accelerated sine wave as described above. The continuum 
curves are also presented. The cycle time T = 60 msec has 
been used for these functions of time. 
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In terms of scene presentation delay, figure 8(b) repre- 
sents the current state of the art in CGI technology. It 
exemplifies the phase problems for frequencies in the 
bandwidth out to about 20 rad/sec. For high frequencies 
the projection operation of the P-E computer amplifies 
signals slightly, but this amplification is minimal, especially 
when compared to the amplification that accompanies 
compensation algorithms. 

The projection operation of the P-E computer is critical 
in minimizing the granularity caused by both asynchronous 
operations and large host computer cycle times. The clock- 
measured projection interval T ^ is as given in both figures 
1 and 6(b). Note, however, that the applicable position and 
velocity signals are modified within the host computer when 
the compensation algorithm is used. The algorithm does not 
require changes within the CGI system. 

The standard projection operation of the P-E computer is 
used in all of the following material. 

Application of the Algorithm 

By using the compensation algorithm, decreases in phase 
error (or time delay) are accompanied by alterations in 
magnitude responses. This is demonstrated below by varying 
the available parameter W which has been defined as the 
cutoff frequency. A satisfactory W selection is dependent 
upon the frequency content of the position and velocity 
data transmitted by the host computer. 

Without compensation, the applicable magnitude and 
phase responses have been given in figures 7(a) and (b) by 
the particular curves labeled “with projection.” The perti- 
nent time response has been given in figure 8(b). These 
curves are for comparison with the following cases where the 
compensation algorithm is used. 

Using the compensation algorithm, figures 9(a) and (b) 
are produced. They give the magnitude and phase variations 
as a function of the parameter W (cutoff frequency). For 
convenience, the continuous system response is also included: 
in these figures. Figures 10(a) and (b) complete the picture 
by giving time responses using two values of the selected 
parameter. From these figures, the tradeoff becomes 
apparent: If the data is band-limited to 10 or 15 rad/sec, 
amplification beyond this limit may be ignored, and the 
algorithm clearly eliminates the objectionable phase error. 
If higher frequencies are known to exist in the data, the 
parameter W must be increased with corresponding degrada- 
tion in the lower frequency region. Although this degrada- 
tion is not rapid, alternate techniques for attenuating the 
high-frequency gain may be required in certain applications. 
Specifically, anti-aliasing or notch filters may be required 
for blade-element rotorcraft models (ref. 12). 

Compensation for time delay is invariably accom- 
panied by magnitude variations. The essential differences 
between this discretely-designed compensation scheme and 


continuously-designed schemes are: l)the concomitant 

reduction to discrete form is avoided, with its resultant 
deviation from the design objective; 2) two distinct fre- 
quencies are used in the design process, with good char- 
acteristics over the associated interval. By defining one 
distinct frequency as d.c. and the other at 10 or 15 rad/sec, 
the region of crossover is generally included (0 < w c < W). 
Hence, the question of design relative to a particular aircraft 
or pilot is avoided. As with analogous compensation tech- 
niques, the “resulting increase in system gain at frequencies 
>w c is not normally a problem, because system amplitude 
ratio and input and disturbance signal power usually decrease 
rapidly at frequencies >w c ” (ref. 9). 

If signal power does indeed decrease rapidly beyond the 
crossover frequency, additional attenuation is not required. 
The single parameter W should of course be small with 
respect to the Nyquist frequency. Since the operational 
bandwidth of the pilot should also be small with respect to 
the Nyquist frequency for successful real-time simulation, 
the optimum placement of W can be at the discretion of the 
simulation analyst. The results have proven to be relatively 
insensitive to this parameter. 

A typical single-axis dynamic check is shown in figure 1 1 , 
where the parameter W = 10 rad/sec. For these particular 
real-time simulation traces the required cycle time was 
62 msec. In figure 11, the positional lead of approximately 
130 msec is observed in the CGI commanded roll angle, as 
predicted. Also, modifications are apparent in the projection 
velocity, as transmitted to the CGI system. 

A novel experiment was conducted to verify both the 
existence of time delays, as quantified herein, and the com- 
pensation algorithm’s ability to eliminate these delays. This 
experiment used sinesoidal commands over a range of fre- 
quencies to drive both an oscilloscope display, and a CGI- 
created horizon. Various host-computer cycle times were also 
used. By use of optical beam -splitting techniques the visual 
signals appeared on the same CGI monitor. Superimposition 
of these signals was possible by manual adjustment of a phase 
component in the oscilloscope drive signal. Researchers 
replicated the standard CGI system’s theoretical delay values 
to within negligible tolerances. Most importantly, they could 
not identify time delays whenever the compensation 
algorithm was used. 

CONCLUSIONS 


A compensation algorithm has been developed that 
virtually eliminates CGI scene presentation delay for the 
frequency bandwidth out to about 15 rad/sec, using a 
worst-case scenario with cycle times in the 60-msec range. 
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If position and orientation signals are known to be band- 
limited by this value, the compensation algorithm is appli- 
cable without further consideration. 

If high-frequency content is known to contaminate the 
CGI drive signals, further considerations may be necessary 
such as anti-aliasing or notch filters. 


In general, the compensation algorithm described here 
should enhance CGI performance, and be especially 
valuable for high-gain tasks such as those required in 
helicopter simulation. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, California 94035, May 14, 1985 
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APPENDIX I 


PERIODIC DELAY 


The periodic component of CGI time delay may be 
written 

T^(m) = T^mM/N - [mM/N\\ 

where the half bracket operation denotes the least integer. In 
this equation the ratio of sample rates 

f/T= 1/(30 T)=M/N 

has been expressed in terms of M and N that are relatively 
prime. The average value for the periodic component of time 
delay is given by 

r /i (av8) “| 2 r <4 ( "° 

m = l 

which may be evaluated using the following relationships. 
First, when the greatest common divisor of the positive inte- 
gers Af and N is unity we have (ref. 1 2) 

N 

£ [mM/N\ =M+ (M- IXN- l)/2 
m = 1 


and second, we have the trivial summation 
N 

J^m=N(N+ l)/2 
m = 1 

These produce the average periodic delay: 

N 

T^Cavg) = (T/N) £ {mM/N - [mMjN\ ( 
m = 1 

= (T/N){(M/N)[N(N + l)/2 - [M{\ +N)+ 1 -iV]} 

= (r/ 2 )(i - 1 /ao 

Also, it may be demonstrated that for m on (1 ,70) the maxi- 
mum value for T^(m) is twice the average value. It is conjec- 
tured that for all m, 

7^ (max) = 7X1 - 1/AO 

which means that the total time delay is never quite as large 
as the sum of cycle time and pure transport delay. 
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Figure 1 Asynchronous data flow. 
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Figure 3.- Scene presentation delays, (a) Cycle time of 49 msec, (b) Cycle time of 50 msec, (c) Cycle time of 51 msec 

(d) Cycle time of 52 msec. 
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Figure 3.- Concluded, (i) Cycle time of 57 msec, (j) Cycle time of 58 msec, (k) Cycle time of 59 msec. (1) Cycle time of 

60 msec. 
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Figure 4 - Continuous and discrete transfer function, (a) Magnitude, (b) Phase. 
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Figure 5.— Transfer function and linear compensation, (a) Magnitude, (b) Phase. 
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Figure 8.— Time history of system without compensation, (a) Without projection, (b) With projection. 
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Figure 10.- Time history of system with compensation, (a) Cutoff frequency of 10 rad/sec. (b) Cutoff frequency of 15 rad/sec 
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Figure 1 1 Typical dynamic check, roll angle and rate with compensated CGI commands. 
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